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Traditionally stellar radiation has been the only heat source considered capable of deter- 
mining global climate on long timescales. Here we show that terrestrial exoplanets orbiting 
low-mass stars may be tidally heated at high enough levels to induce a runaway greenhouse 
for a long enough duration for all the hydrogen to escape. Without hydrogen, the planet no 
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longer has water and cannot support life. We call these planets "Tidal Venuses," and the 
phenomenon a "tidal greenhouse." Tidal effects also circularize the orbit, which decreases 
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tidal heating. Hence, some planets may form with large eccentricity, with its accompanying 
large tidal heating, and lose their water, but eventually settle into nearly circular orbits 
(i.e. with negligible tidal heating) in the habitable zone (HZ). However, these planets are 
not habitable as past tidal heating desiccated them, and hence should not be ranked highly 
for detailed follow-up observations aimed at detecting biosignatures. We simulate the evo- 



lution of hypothetical planetary systems in a quasi-continuous parameter distribution and 
find that we can constrain the history of the system by statistical arguments. Planets or- 
biting stars with masses <0.3 Ms un may be in danger of desiccation via tidal heating. 
We apply these concepts to Gl 667C c, a ~ 4.5 Msarth planet orbiting a 0.3 Ms un star at 
0.12 AU. We find that it probably did not lose its water via tidal heating as orbital stability 
is unlikely for the high eccentricities required for the tidal greenhouse. As the inner edge 
of the HZ is defined by the onset of a runaway or moist greenhouse powered by radiation, 
our results represent a fundamental revision to the HZ for non-circular orbits. In the ap- 
pendices we review a) the moist and runaway greenhouses, b) hydrogen escape, c) stellar 
mass-radius and mass-luminosity relations, d) terrestrial planet mass-radius relations, and 
CN ' e) linear tidal theories. 
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1 Introduction 



All life on Earth requires liquid water. Not surprisingly, o ur search for life on exoplanets 



therefore begins with identifying environ ments that support it (IHartl . Il979t iKasting et al.l . 11993 



Selsis et all 120071 : Ivon Bloh et all 120071 1 . Many other features of a planet are also important, 
such as the mix of gases in the atmosphere and the interior's structure and energy budget, but 
current research has focused primarily on the stability of surface water. On Earth, liquid water 
persists primarily because the Sun's radiation heats the surface to a temperature between water's 
freezing and boiling points. Thus, the concept of a "habitable zone" (HZ) emerged, which is the 
region around a star in which insolation can maintain liquid water on the surface, assuming an 
Earth-like planet. 

Stars have a wide range of luminosities, most of which are considerably lower than the 
Sun's, yielding HZs that are closer in. At the extreme, the luminosity becomes so low that 
a planet orbiting at the stellar surface would be too cold to support liquid water on its sur- 
face. Such a star has no HZ. Th e recently-discovered Y dwarf (which is not a star) WISEP 
J1828+2650 ( jCushing et all 1201 ll ). with an effective temperature below 300 K, is an example of 
such a primary. For warmer stars, HZs may still be cl ose enough in that non-radiative processes 
may i mpac t habi t ability, such as stel l ar flaring (e.g. iLammer et all 120071 ; iKhodachenko et al. 



200 71: iTianl. l2009t ISeerura et al 



2010t ILammer et all 120100 . d ecreased initia l volat i le inventory 



( Ra ymond et all 120071: iLissauerl. I2007Q or tidal effects (e.g. IKasting et al.l. Il993l; iJoshi et al. 



19971 : iJackson et all l2008at ICorreia et all 12008 : Barnes et all l2009al : iHeller et all 1201 If ). As the 

HZ of our Sun is too distant for these phenomena to affect the Earth, we can currently only 
explore their role theoretically, and, consequently, many scientists consider close-in planets less 
favorable candidates for habitability. 

Nevertheless, the last few years have see n renewed interest in the potential habitability of 



planets in orbit about low luminosity objects (jTarter et all l2007t IScalo et all 120071 : iLunine et al. 



20081 : iMonteirol . I201CH ; lAgoll . 1201 It iBolmont et all 1201 if ). This shift occurred because terrestrial- 
sized planets are easiest to detect around low luminosity hosts, due to the larger mass and radius 
of the planet relative to the star. Furthermore, these objects are the most abundant in the 
solar neighborhood. Planets with masses between 1 and 10 Msarth, "super-Earths," have indeed 
been detected in t he last few years around low- mass stars, such as Gl 581 d and Gl 667C c 



Bonfils et al.. 2011 



by radial velocity (Udrv et all 120071: iMavor et all. 120091: IVoet et al.l. 120101 iForveille et all 12011 



Anriada-Escude et all 120121: iDelfosse et all I2012J ) and GJ 1214 b by transit 



( jCharbonneau et all 120091 : ISada et all 120101 : Kundurthy et all l2011t ICarter et all 120111 ) . Sev- 



eral observation al campaigns designed specif i cally t o detect planets around low mass st a rs are 



Bean et all l2010t iRodler et all J201lh . 



now underway (INutzman and Charbonneaul . 120081 : iBoss et all l2009t IZechmeister et all 12009 



The inner edge of the HZ is especially prone to non-radiative phenomena. Traditionally, the 



inner boundary is denned by the stellar dista nce at which in solation is strong enough to remove all 
water in the atmosphere and on the surface (jKastingi . Il988l ) . Two different dehydrating scenarios 
(the "moist" and "runaway" greenhouse) are discussed in detail in §[2] and App. |A] Both require 
water vapor to penetrate a stratosphere, be dissociated (photolyzed) by high energy radiation, 
and culminate in the escape of hydrogen. We call any process that can ultimately lead to total 
water loss a "desiccating greenhouse." Without sufficient hydrogen, water cannot form, and the 
planet will remain uninhabitable forever, unless a major, unlikely event occurs, e.g. an impact 
by a water-rich body that simultaneously delivers water and changes the atmosphere in a way 
that halts the desiccating greenhouse. This definition of the inner edge is conservative because 
no known process can maintain habitability against a desiccating greenhouse. Other processes 
could be equally deleterious for life, but clearly total desiccation will terminate habitability. 

We will discuss four types of terrestrial planets in this study, classified by their water content. 
"Wet" planets are terrestrial exoplanets that have a water content similar to the Earth. "Dry" 
exoplanets have far less, ~ 3 cm deep if condensed and spread globally, but are habitable (see 
Abe et all 120111 ) . "Desic cated" planets have lost all their water th rough a desiccating greenhouse. 
Finally, "water worlds," (JRaymond et all I2004J ; iLeger et al.l . l2004l ) are planets with a much larger 
inventory of water than the modern Earth, e.g. a warm Europa. This investigation explores the 
tra nsition of w et, dry and water worlds into the desiccated state. 



Dold (11964] ) was the first to point out that terrestrial planets in the HZ of low luminosity stars 



can have their spin altered by tidal interaction. In particular , the danger of synchr onous rotation, 
i. e. one hemisphere always facing the star, was emphasized. iKasting et al.l (119931 ) quantified this 
concept and found that planets orbiting within the HZ of stars less than two-thirds the mass 
of the Sun were in danger of synchronization. Although their analysis was limited to Earth-like 
planets on circular orbits, a general belief developed that those planets could not be habitable, 
as one half of the planet would freeze while the other would roast. 

Synchronization is certainly an important consideration when assessing habitability, and 
many investigations have explored its role, but with mixed results. Atmospheric modeling ini- 
tially suggested that cir culation will trans port energy to the unlit side, ameliorating the extreme 
temperature difference (IJoshi et all 119971 ). Some subsequent mod eling has confirmed that syn- 



chronous rotators are likely to have super-rotating atmospheres (jJoshil . 120031 ; iHeng and Vogt 



2011 



Edson et all l201ll ; IShowman and Polvanil . l201ll ) , while others have discounted it (IWordsworth et al. 



201ll ). a nd still others have su ggested that such a state might be beneficial at the outer edge of 
the HZ (jPierrehumbertl . 1201 ll ). Taken together, these investigations suggest that synchronized 
planets should not be dismissed uniformly as uninhabitable. Hence, synchronization is not as 
stringent a constraint as desiccation, and therefore is not an HZ boundary. 

For many years, confusion also existed regarding the term "tidal locking." Many investiga- 
tors assumed it was syn onymous with "syn chronous rotation." If the orbits are non-circular, 



as for many exoplanets (JButler et al.l . I2006J ). then tidally-evolved planets may reach an equi- 



librium state where they rotate faster than synchronous with an "equilibrium" o r "pseudo 



synchronous" period. This aspect of tidal theory has been known for decades (e.g. iGoldreichl . 



19661 ; iGreenberg and Weidenschillind. I1984T). but has only rec e ntly been poin t ed ou t for the case 



of exoplanets f lBarnes et all 120081 : iFerraz-Mello et all 120081 : ICorreia et all 120 08). 



some exoplanets, such as Gl 581 d with an ecce ntricity of 0-3 8 (IMavor et al. 
"tidally locked" but rotat e about twice per orbit (jBarnes et all l2008t iHeller et al 
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201lh . For 



more on this point, consult IMurray and Dermottl (J1999J), chap. 5.2, and for a counterargument see 
Makarov and Efroimskyl ( 120121 ) . In this paper, we use "tidally locked" to mean a planet rotating 
at the equilibrium period as determined by its eccentricity and obliquity, see Eqs. ( 123]) . (I3~TJ) -( 132]) . 
In summary, a synchronously rotating planet is tidally locked (yet, it could still have non-zero 
eccentricity and obliquity), but a tidally locked planet is not necessarily rotating synchronously. 
Even if an orbit is currently circular, tides may not have driven the rotation to synchronous. 
If the orbit began with large eccentricity, tides will tend to damp it to zero and we may expect 
it to be rotating synchronously. However, the planet could pass through one or more "spin-orbit 
resonanc es," where the planet's rotational frequency is commensurate with its orbital frequency 



(see e.g. 



Rodriguez et all 120121 ). For example, Mercury rotates three times for every two times 
it orbits the sun, a 3:2 spin-orbit resonance. Spin-orbit resonances require an inhomogeneous 
mass distribution, which is likely for tidally-deformed exoplanets, but cannot be measured for 
the foreseeable future. A planet caught in a spin-orbit resonance may remain in that state even if 
circularized, as a resonance is a strong dynamical process. For any particular exoplanet, capture 
and retention into a spin-orbit resonance will be very difficult to constrain observationally, so all 
reasonable options should be considered. Therefore, synchronous rotation is unlikely for planets 
with large eccentricities (Mercury's is 0.2), and not e yen guaranteed for a circula r orbit. For 
more on spin-orbit resonances, the reader is referred to IMurray and Dermottl ( 119991 ). chap. 5.4. 
As tidal locking of the planetary rotation is not an absolute constraint on habitability, we 
turn to tidal heating as the other tidal phenomenon most likely to affect planetary habitability. 
As a planet moves from periastron, its closest approach to the star, to apoastron, the furthest 
point, and back again, the gravitational force changes, being inversely proportional to distance 
squared. This difference creates an oscillating strain on the planet that causes its shape to vary 
periodically. The rigidity of the planet resists the changes in shape, and friction generates heat. 
This energy production is called tidal heating. 



Tidal heating is responsible for the volcan ism on Io (IStrom et all Il979t lLaver et all 120071 ). 



which was predicted, using tidal theory, by IPeale et all (119791 ) . Io is a small body orbiting 



Jupiter with an eccentricity of 0.0041, which is maintained by the gravitational perturbations of 
its fellow Galilean moons, that shows global volcanism which resurfaces the p l anet on a timescale 



of 100 - IO 5 years (IJohnson et all 11979 : iBlanev et all 1199a McEwen et all 120041 ) . The masses 



of Jupiter and Io are orders of magnitude smaller than a star and terrestrial exoplanet, and 
thus the latter have a much larger reservoir of orbital and rotational energy available for tidal 



heating. Moreove r , some exop l anets have been found wi th orbital eccentricities larger than 0.9 



(INaef et all l2001t Ijones et all 120061 ; iTamuz et all . 20081). Thus, the tid al heating of terrestrial 



exoplanets may be much more effective than on Io ( Jackson et all . l2008d Jal; iBarnes et all l2009al . 



20101 ; iHeller et all 1201 ll ). This expectation led to the proposition that terrestrial exoplanets with 
surface heat flux es as large or larger than Io's should be classified as "Super- Ios", rather than 



"Super-Earths" (IBarnes et all |2009bJ) . Numerou s Super-Io candidat es exist, su ch as CoRoT-7 b 
dLeger et all 120091: IBarnes et allboiol), G1581 e (TMavor et alll2009h . 55 C nc e dMcArthur et al " 



2004J : iDawson and Fabrvckvl . l2010l : IWinn et all |201l|), and Kepler- 10 b featalha et all l201lh . 

but none is in the HZ. 



Jackson et al . (2008ai) a ndlBarnes et al.l (I2009at) considered Io's heat flux, ~ 2 W m 2 (IVeeder et al. 



1994 : ISpencer et all I2OOO1 : iMcEwen et all . 120041 ). to be an upper limit for habitability, arguing 
that Io-like surfaces are dangerous for habitability. However, 2 W m -2 may not be sufficient to 
sterilize a planet, and so should not be considered a hard limit to habitability. For example, the 
heat flow in inhabited hy drothermal vent sys tems on Earth, suc h as the Endeavour s egment of 



the Juan de Fuca Ridge (IHolden et al.l . Il998l ). is ~ 30 W m 2 (IFontaine et al.l . l201ll ). Thus, a 



water world with tens of W m -2 of energy output could support life. While Io-like volcanism 
is clearly an issue for habitability, it may not always lead to sterilization. However, if the tidal 
heating can maintain a desiccating greenhouse long enough for all the water to be lost, then the 
planet becomes uninhabitable, and is highly unlikely to ever regain habitability. 

Calculations of the inner edge of the HZ have traditionally assumed the primary energy so urce 
at the surface is stellar radiation, as is the case for the Earth. Following IBarnes et al.l (l2009al ). we 
call that type of HZ an "insolation HZ" (IHZ), as starlight is the only energy source considered. In 
this investigation, we identify the amount of tidal heating that triggers a desiccating greenhouse, 
as well as the combinations of physical and orbital parameters for which tidal heating could yield 
it. We find such a state is predicted by current models, and dub such a world a "Tidal Venus," 
and the phenomenon that produces it a "tidal greenhouse." The tidal greenhouse would probably 
have the same effect on habitability as one caused by irradiation, and hence should be considered 
as hard a limit to habitability as a traditional, insolation-driven desiccating greenhouse. 

In this study, we define the limits of Tidal Venuses in terms of stellar and planetary mass, 
M* and M p , respectively, the orbital semi-major axis a, orbital eccentricity e, planetary radius 
R p , planetary obliquity ip p , and planetary spin frequency uj p . Some of these quantities are easily 
observed by current technology, others require special geometries and the next generation of 
space telescopes. Therefore, for the next decade, not all newly-found planets in an IHZ will have 
well-constrained tidal heating. This study provides a framework for identifying the range of tidal 
heating on terrestrial planets orbiting in the IHZ of low luminosity stars. As we see below, the 
story is complicated and involves a large parameter space, but is tractable. 

We first review the surface conditions that lead to a greenhouse state in § [2], including an 
extended discussion of the moist and runaway greenhouses (App. IA~|) . a short review of hydrogen 



escape driven by high energy radiation (App. [B]), the relationships between mass, radius and 
luminosity for low-mass, hydrogen burning stars, "M dwarfs" or "late-type stars" in the parlance 
of astrophysics, the extent of the IHZ (App. [C]), and mass-radius relationships for terrestrial 
exoplanets (App. [D]). We then briefly describe tidal heating in § [3] with details in App. [El Next 
we show that tidal heating alone can produce surface conditions capable of triggering a runaway 
greenhouse on planets orbiting M dwarfs (§ H]). Then we explore how past tidal heating may 
preclude habitability of planets found in the IHZ regardless of their current tidal heating (§ E]). 
In § [6] we then consider the Gl 667C system which contains two potentially habitable planets 
and find that tidal heating is unlikely to have sterilized either. In §[7] we discuss the results, and 
finally, in § [8] we draw our conclusions. 

2 The Inner Edge of the Habitable Zone 

Planets with surface water will always have water vapor in their atmospheres. The amount of 
water vapor present in equilibrium with a liquid surface is described by the saturation vapor 
pressure p sat , which is the pressure exerted by water vapor in thermodynamic equilibrium with 
standing surface water. The value of p sat depends exponentially on temperature: a warmer planet 
will have much more water vapor in its atmosphere. Water vapor is a greenhouse gas (see below), 
so a warmer planet will have a stronger greenhouse effect, enhancing warming, a positive (but 
not necessarily runaway) feedback. 

Atmospheric gases on Earth — and those that we expect on habitable Earth-like planets — are 
mostly transparent to optical wavelengths, so starlight is able to heat the surface. The surface 
in turn heats the air in contact with it, which rises, expanding and cooling adiabatically. As the 
air cools, water vapor condenses, releasing latent heat and so slowing the cooling. The net rate 
of temperature decrease with height is called the "moist adiabatic lapse rate," setting the mean 
temperature-pressure (T-p) structure of Earth's troposphere (the lower region of the atmosphere 
affected by convection, and bounded at the bottom and top by the planetary surface and the 
tropopause, the altitude at which temperature stops dropping with height). 

Atmospheric gases which absorb radiation at similar wavelengths to the radiative emission of 
the planet (thermal infrared for Earth and any habitable planet) are termed "greenhouse gases" ; 
these cause a "greenhouse effect". The most important greenhouse gases for Earth are water 
vapor and carbon dioxide. They absorb radiation emitted by the surface, and then re-emit, 
both downward toward the surface and upward into space. Energy radiated towards the surface 
(commonly called back radiation) heats the surface. Most importantly, as the atmosphere is 
cooler than the surface, the amount of radiation emitted to space from the atmosphere is less 
than the amount of energy emitted by the surface. Thus, the presence of a greenhouse atmosphere 
means the surface temperature is warmer than the effective temperature of the planet. 



Earth's greenhouse effect keeps the planet warm enough to be habitable; without it Earth's 
surface would be at the planet's effective temperature, a barren 255 K. However, stronger and 
more water-vapor-rich greenhouse atmospheres can render the planet uninhabitable through high 
temperature sterilization and/or desiccation (loss of the ocean), and thus define the inner edge of 
the habitable zone. There are two physically distinct situations that lead to loss of habitability, 
the "runaway greenhouse" and "moist greenhouse" . 



The runaway greenhouse was recently reviewed by iGoldblatt and Watson! (120121 ) . so we pro- 
vide a summary description only. As the planet warms, the amount of water vapor in the 
atmosphere increases such that water becomes a major constituent of the atmosphere and ul- 
timately the dominant one. Consequently, the moist adiabatic lapse rate tends toward the 
saturation vapor pressure curve for water and the T-p structure of the atmosphere becomes 
fixed. Concurrently, the atmosphere becomes optically thick in the thermal infrared, such that 
only the upper troposphere can emit to space. As the T-p structure is fixed, the emitted ra- 
diation is also fixed, imposing a limit on the outgoing radiation (F cr u) from the troposphere. 
Values for F rrit from the literature are typically 285 Wm~ 2 to 310 Wm~ 2 for a 1 MKarth plane t 
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flPollackl . 119711 : IWatson et all I1984J ; lAbe and Matsuil . ll988l:lKastind. Il988l: Pierrehumbcrt 



2010). 



The p h ysics of this limit is desc ribed in more d etail in iSimpsonl (119271 ) and iNakajima et al. 
(119921 ). (JKomabayashil (119671 ) and llngersolll (119691 ) describe a stratospheric limit at a higher flux, 
385 Win -2 , which is never reache d in practi ce.) If the amount of e nergy supplied to the atmo- 



sphere by the Sun (ISimpsonl . Il927l ). impacts (jAbe and Matsuil . Il988l ). or tidal heating (this work) 
was to exceed F crit then the atmosphere would not be able to maintain radiation balance and 
runaway heating of the surface would ensue, causing evaporation of the entire ocean. Radiation 
balance would be regained when either a) the surface temperature reaches ~ 1400 K at which 
point enough radiation is emitted in the near infrared where water vapor is not a good absorber, 
or b) if all the water vapor is lost from the atmosphere. 

One likely water loss mechanism is hydrogen escape to space. Today, little water vapor 
reaches the upper atmosphere because the tropopause acts as a "cold trap" : almost all water has 
condensed lower in the atmosphere, so the tropopause water vapor mixing ratio (the abundance 
of water vapor in the atmosphere, expressed as the ratio of the mass of water vapor to the mass 
of dry air) is low. Water vapor transport above here is generally diffusive, leading to a constant 
water vapor mixing ratio in the upper atmosphere. (More precisely, in Earth's atmosphere, water 
vapor increases with altitude in the stratosphere due to methane oxidation and decreases with 
altitude in the next higher atmospheric layer, the mesosphere, due to photolysis; however, the 
total hydrogen mixing ratio is conserved, and the hydrogen escape rate depends on this value.) 
Thus, water is not a strong source of escaping hydrogen and the overall escape rate is low. 
However, in the runaway greenhouse the atmosphere becomes predomina ntly water, so no such 
const raint applies and hydrogen may escape hydrodynamically to space (JKasting and Pollack! . 
19831 ). The rate will depend on the amount of extreme ultraviolet (XUV, 1 - 1200A) radiation 



absorbed in the hig hest part of the atm osphere, so this process is also limited by the available 
stellar XUV energy (IWatson et al.L Il98lh . A Venus-size planet at the inner edge of the habitable 



198ll ); we refer 



zone, could have lost an ocean the size of Earth's in ~ 10 years (IWatson et al 

to this as the desiccation time, tdes (discussed further below). Deuterium would be retained 

preferentially over ordinary hydrogen during hydrogen escape. Enrichment o f D/H on Venus 



19821 : Ide Bergh . 



impli es that Venus lost a substantial amount of water to space (jDonahue et al 
19931 ). likely after experiencing a runaway greenhouse. 

We adopt td es = 10 8 years based both on the work on Venus, as well as calculations presented 
in App. [B] for a terrestrial planet orbiting a very low mass star. Suc h stars are very a ctive 



and emit strongly in the near UV, especially when th ey are young (e . g. I West et al.l . 120081 ). but 



the total flux is still less than the present-day Sun ( JFleming et al.l . Il993l ) and the XUV flux 



(which drives hydrodynamic escape) is not well known. Thus, we expect a range of td es to exist 
especially since the orbit evolves with time due to tides, but choose this value for simplicity. 
The model described in App. [B] can be applied to individual cases when observations permit 
realistic modeling. We stress that there are many unknowns which may affect the actual value 
of tdes- The ma sses of exoplanet oceans are unknown. Earth's mantle contains a few ocean 



masses of water (JBell and Rossmanl . Il992l ; iMurakami et al.l . |2002| ; iMarty and Yokochiml . 120061 ). 



so it is possible that the surface of a planet could acquire a new ocean via outgassing after 
desiccation, and become habitable again. Thus, we cannot conclude with certainty that a planet 
with a desiccated atmosphere will never be habitable in the future. However we also expect rapid 
mantle overturn during strong tidal heating, and hence the entire water inventory may actually 
be lost. We therefore assume that the tidal greenhouse is highly likely to permanently sterilize 
a planet. 



The moist greenhouse ( IKastingl . Il988l ) describes a warm atmosphere, in which the whole tro- 
posphere is assumed to be water vapor saturated and is underlain by a liquid ocean. As the 
surface temperature increases, the tropopause is pushed higher (to lower pressure) given the rea- 
sonable assumption of a constant tropopause temperature. While the saturation vapor pressure 
at the tropopause is, by this assumption, constant, the saturation mixing ratio of water at the 
tropopause (p sa t/ptro P ) increases as the tropopause moves higher, where pt Top is the atmospheric 
pressure at the tropopause. The cold trap is then no longer effective, and substantial water vapor 
penetrates the stratosphere, and water-derived hydrogen escape can be effective. The planet may 
thus gradually desiccate while in a hot, but stable, climate. This process could be driven by high 
greenhouse gas inventories rather than external heating. 

A key distinction between the runaway and moist greenhouses is that the former happens 
when a known flux of energy is supplied to the planet, whereas the latter depends most strongly 
on surface temperature. The runaway greenhouse occurs because water vapor is a greenhouse 
gas and can therefore trap heat near the surface, and hence is only a function of the absorptive 
properties of water and the energy flux. The moist greenhouse occurs when the temperature at 
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the tropopause is large enough that water does not condense and is therefore able to escape to 
the stratosphere where it can be photolyzed. Increasing the planet's inventory of non-condensible 
greenhouse gases has a small effect on the runaway greenhouse limit, but can drive the moist 
greenhouse since the relative amount of H2O is lower. Thus, the runaway greenhouse is a more 
conservative choice to demarcate the inner edge of the IHZ in the sense that it depends solely 
on water, whereas knowledge of the atmospheric composition may not be available to assess the 
likelihood of a moist greenhouse. 

In Fig. [T]we show the moist and runaway greenhouse for wet planets in orbit around M dwarfs. 
See Apps. |Aj [C] and |D] for a discussion of IHZ limits, stellar mass-radius and mass-luminosity 
relationships, and terrestrial mass-radi us relationships. In Fig. [TJ the grey re gions are the limits 
of the moist greenhouse as presented in lSelsis et al.l (120071 ) . At the inner edge, ISelsis et al.l (120071 ) 
find that 300 W m -2 triggers the moist greenhouse on a 1 M Earth planet, and, from left to right, 
the limits assumed 100%, 50% and 0% cloud cover (for a 0.25 Ms un star these limits correspond 
to bo nd albedos of 0-75, 0.4 9 and 0.23, respectively). The solid curve is the runaway greenhouse 
limit (IPierrehumbertl . I2OIOI ) for a 30 MEarth planet and dashed for a 0.3 MEarth planet, both with 
an albedo of 0.49 (compare to the medium grey). As expected (see App. |AJ the smaller planet's 
inner edge lies at larger semi-major axis than both 1 ME art h moist greenhouse limit, which in 
turn lies farther out than the 30 MEarth runaway greenhouse limit. 

We have presented the classical description of the desiccation at the inner edge of the habitable 
zone above. Various complications are worthy of note. Water on a dry planet will get trapped at 
the poles, making a moist or runa way greenhouse h arder to achieve and meaning that the inner 
edge of the HZ is nearer the star (jAbe et all l201ll ). We do not include that limit in Fig. [CJbut 
we return to it in § SJ Fig. [T] assumes t he planetary orbit is circular. However, for the many 



exoplanets with (highly) eccentric orbits (IButler et all 120061). the total irradiation over an orbit 



-1/4 



(Barnes et all 120081 ). 



determines the annual-averaged surface temperatu re ([Williams and Pollardl . |2002| ) and pushes 
the IHZ boundaries out by a factor of (1 — e 2 



3 Tidal Heating 



Tidal heating is respo nsible for the volcanic activity on Io (IPeale et all 11979). an d probably the 



geysers of Enceladus (IHansen et al.l . 120061 ; iPorco et all 120061 ; iHurford et all 120071 ). Tidal theory 
has a long and established body of work, however tidal processes remain poorly understood. The 
difficulty lies in the complexity of the energy dissipation processes and the very long timescales 
associated with tidal evolution. The Earth-Moon system is the most accurately studied, with 
lunar laser ra nging; providing a pr ecise measurement of the recession of the Moon due to tides of 
~ 38 mm/yr toickev et al. . 1994), as well as direct measurements of the locations of dissipation 
in the ocean (JEgbert and Ravi . l2000h . However, using the currently estimated tidal dissipation 
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Figure 1: Comparison of HZ boundaries for different planeta ry masses and de siccating green- 
houses. The shaded regions represent the IHZ boundaries from lSelsis et al.l (120071 ): Dark gray as- 
sumes no cloud coverage, medium gray 50%, and light gray 1 00%, and assum i ng 1 M Earth planet. 
The black curves represent the runaway greenhouse limit from iPierrehumbertl ( 120111 ). with a plan- 
etary albedo of 0.49. The solid curve is for a 30 M Earth planet, dashed for a 0.3 Mp, ai . t h planet. Fo r 
these calculations we used the stellar mass-radius relationship from iBayless and Oroszl ( 120061 ) , 
the mass-luminos ity relationship fro m iReid and Hawleyl (120001 ) , and the terrestrial mass-radius 
relationship from lSotin et al.l ( 120071 ). 
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parameters to extrap olate backwards in time predicts that the Moon was at the Earth's surface 
about ~ 2 Gyr ago ( JMacDonaldl . Il964l ). which contradicts the standard impact model for the 
lunar origi n (note that numerous issu es remain such as the origin of the Earth and Moon 's 
obliquities (ITouma and Wisdom! . 119941 ) and the perturbations from other planets (jCukj, 120071 )). 
The deceleration of Io's orbit al velocity has also been tentatively detected a nd seems broadly 



consistent with tidal theory (jAksnes and Franklin! l200ll ; lLainey et al.l . 120091 ). However, these 



data lie in the same regime, that of nearly-circular orbits. Exoplanets have been found with 
extremely eccentric orbits, so we can only extrapolate from our Solar System cautiously. The 
details of tidal theory are complicated, and hence we relegate the discussion to App. [Ej 

Two end-member models of tides theory have been applied to exoplanets and the bodies of 
our Solar System: One assumes the energy dissipation gives a constan t phase lag i n the periodic 
distortion (CPL), and the other assumes a constant time lag (CTL) ( JGreenbergl . 120091 ) . While 
simple and linear, such models are probably commensurate with the dearth of information we 
have about exoplanet interior processes. More complicated models have been constructed, and 
they reproduce t he above mode l s for c ertain choices of internal composition, structure, and energy 
transport (e.g. iHenning et al.l . 120091 ). Thus, the CPL and CTL models can provide important 



and accurate insight into the tidal evolution of exoplanets. 

These models con verge at e = 0, and have been shown to be nearly identical for e <0.2 
( JLeconte et al.l . 120101 ) . and when using Eq. (|33|) . However, for e >0.3, they diverge significantly. 
We urge caution when interpreting results in the upcoming sections which allow e to be as large 
as 0.8. We include this range primarily for illustrative purposes, and as a baseline for any future 
work which may include non-linear effects. 

The tidal heating of a body is provided in Eqs. ( 12 2 p and (130 p for the CPL and CTL models, 
respectively. Averaging the heating rate of the entire planetary surface gives the surface energy 
flux due to tides, Fade- In those equations, the strength of the tidal effects is parametrized as 
a "tidal quality factor" Q (CPL) or "tidal time lag" r (CTL), which are notorious ly difficult to 
measure or esti mate from first p r inciples. The Earth's curren t valu es are Q = 12 ( lYoderl . Il995[ ) 
and r = 638 s (jLambeckl . Il977l ; iNeron de Surgy and Laskarl . Il997l ). respectively. However, as 
noted above, these values predict too short a lifetime of the Moon. This discrepancy has led 
to the notion that the Earth's response to lunar tides has evolved with time, possibly due to 
changing size, shape and seafloor topography of the oceans, affecting the ocean currents response 
to the tidal potential. Measurements of the dry bodies in our Solar System have found that their 
Qs tend to cluster around 100, see App. [Ej 

Recent s atellite observations of the Earth have revealed the locations of tidal dissipation in 
our oceans (jEgbert and Ravi . |2000[ ). Tides force water through shallow seas and straits causing 
energy dissipation. In the open ocean, tidal dissipation is probably a non-linear process in 
which currents are disrupted by seafloor topography. The presence of the ocean provides more 
opportunity for tidal dissipation than on dry planets, and we assume that most dissipation on 
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habitable exoplanets will also occur in their oceans. Therefore, for the following calculations we 
assume modern Earth-like planets with Q = 10 — 100 or r = 64 — 640 s. 



4 Tidal Venuses 

We computed tidal heating rates for a range of planetary, stellar and orbital parameters and 
find that tidal heating can be strong enough on some planets to trigger a runaway greenhouse. 
In Fig. [2] we show the configurations that predict this state aroun d 4 different hypo thetical M 
dwarfs. The IHZ boundaries are the moist greenhouse limits from ISelsis et al.l (120071 ) and with 
the same format as Fig. [TJ The colored curves mark where F ti( i e = F crit . Red curves assume the 
CTL model, blue the CPL with discrete rotation states, see App. [EJl. Solid curves assume the 



Pierrehumbertl ( 1201 ll ) runaway greenhouse model, and dotted the dry world model of I Abe et al. 
( 120111 ). c.f. Fig. [6] in App. |A] For these latter worlds, we choose Q = 1 00 as they do not have 
oceans. Thick lines are for a 10 MEarth planet, and thin for 1 MEarth (jAbe et al.l (120 lit ) only 
considered a 1 M Earth planet). Note that our choice for the relationship between Q and r 
influences which m odel predicts more heating. If we had chosen the same relationship as in 
Heller et al.l (120 111 ), we would have found the CTL model predicts more heating than CPL, c.f. 
their Fig. 5. 

For the lowest mass M dwarfs, the HZ is significantly reduced due to tidal heating. For 
masses larger than 0.25 Ms un , a tidal greenhouse in the IHZ is only possible on large mass 
planets with very large eccentricities. Note that our model predicts a tidal greenhouse at low 
eccentricities where tidal theory is most likely valid. Figure |2] shows that a planet may have a 
climate catastrophe due to tide-driven overheating even if it is far enough from the star that 
stellar radiative heating alone would not preclude habitability. 

The difference between curves representing equal mass planets is due to the frequency depen- 
dence of the CTL model. For the Earth, the frequency is the mean motion of the lunar orbit, but 
a planet at a = 0.035 AU, i.e. the middle of the HZ, orbits in about 1 week. One could formally 
adjust r so that it is equivalent to a Q of 10 ((see e.g. iMatsumura et al.l . l201dl ; iHeller et al.l . l201lh 
and Eq. [33]) and then the curves would lie in a similar location. The unknown tidal response 
of terrestrial bodies and the absence of an unambiguous translation from the CPL to the CTL 
model motivates our choice of adopting present-Earth values for all planets at all frequencies. 

Regions to the left of the colored curves in Fig. [2] can produce Tidal Venuses if the planets 
remain there longer than td es - In Fig. El we show the evolution of two example systems consisting 
of a 10 MEarth planet orbiting a 0.1Ms un star, using both the CPL and CTL models. For the 
CPL case, the planet begins with semi-major axis a = 0.04 AU, and e = 0.3. The planet is 
assumed to be spin-locked and with zero obliquity. This orbit is toward the outer edge of the 
IHZ and with the typical eccentricity of known exoplanets. For the CTL case, everything is the 
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Figure 2: Parameter space of the Tidal Venus. The top left panel is a 0.1 Ms un star, top right 
0.15 Msnm bottom left 0.2 Ms un , and bottom right 0.25 Ms un - The grayscale represents the 



Selsis et al.l (120071 ) IHZ boundaries: From lightest to darkest gray, the cloud coverage is 100%, 
50% and 0%, respectively. The colored curves mark where F ti ^. = F r p.t- R ed curves assume the 



CTL model, blue the CPL. Solid curves ass ume the (|Pierrehumbertl . 1201 ll ) runaway greenhouse 



20111 ). Thick lines are for a 10 M 



model, and dotted the dry world model of ( lAbe et al 
planet, and thin for 1 Msarth- Tidal Venuses lie to the left of these curves. 
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Earth 



same, except the initial a is 0.035 AU, as the Tidal Venus region lies closer to the star. The top 
panel in Fig. [3] shows the evolution of a, the next panels down show e, then insolation F inso i, 
then tidal heat flux Fude, and finally the sum of the insolation and tidal heat flux, F to tai- The 
tidal heat flux generated in the planet in the CPL case has decreased to F crit = 345 W mr 2 at 
200 Myr, and the sum of tidal heating and insolation reaches F cr n at 275 Myr. In the CTL 
case, F tide = F crit at 130 Myr, and F tota[ = F crit at 275 Myr. We conclude that planets orbiting 
low-mass stars on eccentric orbits may experience tidal greenhouse conditions long enough to 
become uninhabitable. 

The CPL model predicts a will increase through angular momentum transfer with the star, 
for which we assume a r otation period of 30 days, a tidal Q of 10 6 , a radius determined by the 



Reid and Hawleyl ( 120001 ) relation, and both the radius of gyration and Love number of degree 
2 are 0.5. However, this orbital expansion is a result of the discrete nature of the CPL model, 
which only includes 4 "tidal waves," see App. lEjl. In Fig. [3] the planet initially spins 3 times per 
2 orbits, and hence e^i = 0, see Eq. (J2"T|) . eliminating one of the terms in da/dt. Note that once 
the eccentricity drops below yl/19, the spin rate becomes synchronous (see App.[El3), and then 
£o,i — but e^i = 1, and then da/dt < 0. For this example, that transition occurs at 375 Myr, 
and is not shown in Fig. [3j This point illustrates the complexity inherent in this commonly-used 
tidal model. 

5 Constraining Observed Planets 

The previous section demonstrated that planets may become desiccated by tidal heating, but in 
many cases we will be more interested in the possibility that such a condition developed on a 
planet that we have discovered. As tides tend to circularize orbits, we may find a planet in the 
IHZ with low eccentricity that experienced a tidal greenhouse early on and is hence currently 
uninhabitable. In this section we describe how to evaluate a known planet's probability for 
habitability based on past tidal heating. 

In order to model the tidal heating history of an exoplanet, we require knowledge of M p , R p , 
a, e, M#, R* and age. Exoplanets are predominantly discovered by RV and transit studies that 
can provide a and e, while the others can be modeled from the stellar spectrum. If some of these 
values are unknown, we can often use scaling relations to estimate them, see App. ICHD1 From 
this information, one can estimate the tidal heating history of a planet, using the tidal evolution 
equations (App. |Ej), and hence estimate the probability that the planet has lost its water. 

As an example, consider the hypothetical situation in which a terrestrial-scale planet has been 
discovered in orbit about an M dwarf. Further, the planet lies near the inner edge of the IHZ, 
and has a low eccentricity. To evaluate past tidal heating, we created 4 systems with different 
stellar and planetary properties, applied reasonable uncertainties to each observable properties, 
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Figure 3: Evolution of a lOMEarth planet orbiting a O.lMgun star with an initial orbit of a = 
0.04 AU, e = 0.3. Top: Semi-major axis evolution. Top middle: Eccentricity evolution. Middle: 
Insolation evolution. Bottom Middle: Tidal heat flux evolution. Bottom: Total surface heat flux 
(insolation + tidal). 
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chose non-observable properties via appropriate scaling laws, and modeled the system's history. 
In particular, we assume that the planet's mass lies between 1 and 5 M Earth , has an eccentricity 
less than 0.1, a tidal Q p (r p ) in the range 10 - 100 (640 - 64 s), and is tidally-locked. We assume 
the star's age lies between 2 and 8 Gyr, and Q* (t*) in the range 10 5 - 10 10 (1 - 10~ 4 s) and 
tdes = 10 8 years. We used the same scaling relations as in § [U For the 0.1 Ms un case, a was 
chosen in the range [0.029,0.031] AU, for 0.15 M Sun it was [0.055,0.065] AU, for 0.2 M Sun it was 
[0.085,0.095] AU, and for 0.25 M Sun it was [0.095,0.105] AU. For each parameter the uncertainty 
distribution was uniform in the quoted range, except for e which was chosen uniformly in the 
range —5 < log 10 (e) < — 1. We then randomly determined the system parameters in these ranges 
and integrated the tidal history backward in time, using the models presented in App. [El for the 
randomly chosen age of the system. We ignore the possibility of spin-orbit resonance capture, 
which would dramatically alter the history. For each stellar mass and tidal model, we simulated 
30,000 possible configurations. 

In Fig. H] we show our results graphically. As we are considering a range of masses, we choose 
to represent the runaway greenhouse flux with a 2.5 Msarth planet, i.e. the solid curves show 
where F ti d e = F crit = 309 W m~ 2 . The three contours show the probability density for the 
planet's location td es = 10 8 years after the system's formation. The contours show where the 
probability has dropped by 50%, 90% and 99% of the peak value, solid contours are for the 
CPL model, dashed for the CTL. If the probability contours intersect, or even come close to, the 
colored F U d e = F C rit curves, then the planet may be a "Habitable Zone Venus," a planet that 
appears habitable by the IHZ metric, but is probably more Venus-like than Earth-like. 

The shapes of the regions in Fig. H] are due to our parameter choices and the assumptions 
implicit in the two tidal models. The CPL model does not include eccentricity terms to as high 
an order as the CTL model, and hence the CTL model predicts more evolution at large e. In 
the top left panel, this effect is seen clearly, as the most likely orbits at 100 Myr in the CPL 
model are at lower semi-major axis than the CTL model. We stress that this example is purely 
hypothetical, and the actual shapes of the contours could be very different for actual systems. 

For the 0.1 Ms un case (top left), the planet may have spent enough time in the tidal greenhouse 
to be uninhabitable. For larger mass stars, however, the danger of tidal desiccation is smaller. 
In the 0.15 Ms un case, the peak in the CPL probability density at a = 0.08, e = 0.45 represents 
1% of our simulations. Although the contours do not cross the CPL runaway greenhouse (blue) 
curve, they do come close, and hence there is a small chance that our putative candidate is a 
super- Venus, especially if we allow for absorption of stellar radiation. 

However, for M* > 0.15 Ms U n, planets with low eccentricity probably have always had low 
eccentricity, i.e. the evolution was negligible. This sharp contrast between 0.1 and 0.2 Ms un 
occurs because of the steep dependence of tidal heating on a. At larger stellar masses, the 
circular IHZ has been pushed away from the reach of fatal tidal heating. 

This experiment is purely illustrative with arbitrary parameters and uncertainty distributions. 
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Figure 4: Orbits of hypothetical planets around M dwarfs after t<i es . The top left panel is a 
0.1 M Sun star, top right 0.15 M RllT1 botto m left 0.2 M Sun , and bottom right 0.25 M Sun . The 
grayscale represents the (ISelsis et all 120071 ) IHZ boundaries with the same format as Fig. [2J For 
reference, the red and blue curves show the tidal greenhouse limit for a 2.5 M Earth planet with 
Q p = 10 (CPL model, blue curve) or t p = 640 s (CTL model, red curve). Contours denote levels 
of constant probability density (for densities 50%, 90% and 99% of the peak value) for the initial 
orbit of the planet: solid corresponds to the CPL model (compare to blue curve), dashed to CTL 
(compare to red curve). In the bottom two panels there has been negligible orbital evolution for 
either tidal model. 
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If a terrestrial planet was discovered in the IHZ of a very low mass star, then this methodology 
could be performed in order to characterize its potential to support life. One could derive a value 
of tdes tailored to the primary, and determine the probability that the planet spent more than 
that in a tidal greenhouse state and hence is uninhabitable. In principle, one should also allow 
a range of albedos and include radiation to determine the amount of time that the total surface 
energy flux, F tot = F tide + F insoh is larger than F crit , but we leave such an analysis for future 
work. 



6 Application to Gl 667C 



Two or more planets orbit the 0.3 M Sun star Gl 667C ( jBonfils et all l2011t lAnglada-Escude et al. 



20121 ; iDelfosse et all 120121 ). Planets c and d appear to lie in the IHZ, while a third planet, b, lies 
interior. Planet c is at least 4.5 Msarth and orbits near the inner edge of the IHZ. Planet d, which 
is weakly detected in both studies, lies near the 50% cloud cover outer boundary. At 0.23 AU, 
planet d is too far from the star to be s ubjected to strong tid al heating, at least if it is tidally 
locked. lAnglada-Escude et al.l (120121 ) and IDelfosse et al.l ( 120121 ) propose different solutions to the 
system with the former setting c's eccentricity to 0, but stating that it is only constrained to be 
< 0.27, while the latter assign its eccentricity to be 0.34 ± 0.1 . The minimum mass estimates 



are a lmost identical at 4. 25 Mpwtv, (lAnglada-Escude et al.l 120121 ) and 4.5 Msarth ( IDelfosse et al. . 
2012). Here, we use the lAnglada-Escude et al.l (120121 ) solution, but note that using data from 



( IDelfosse et al.l . 120121 ) does not change our results. 



In Fig. El we show the system in the same format as Fig. |5J As this planet was detected via 
radial velocity data, its true mass is unknown. We consider two possibilities here, an edge-on 
geometry in which its actual mass is 4.5 Msarth and an inclined case in which its actual mass 
is doubled, 9 MEarth- In Fig. [5j the thin lines correspond to the minimum mass, thick to twice- 
minimum. The vertical extent of the line corresponds to the uncertainty in eccentricity. The 
orbit of c is marked by the vertical black line at a = 0.123 AU, and d at 0.23 AU. 

The CTL model (red curves) barely intersect the 100% cloud cover IHZ at large e. The CPL 
model penetrates the IHZ more significantly, but still does not reach the orbit of c. Thus, planet 
c is not currently experiencing a tidal greenhouse, assuming no obliquity and pseudo-synchronous 
rotation. We also see that d is completely safe from the tidal greenhouse. 

But could planet c have been in the tidal greenhouse and be desiccated today? The answer is 
almost assuredly no. In this case, we can appeal to the orbital architecture, rather than running 
a suite of Monte Carlo simulations as in § [5j The dashed, black line of Fig. [5] showed where 
the orbit of c crosses the orbit of b. Such a configuration is as close as the two planets could 
possibly be and remain stable, assum i ng no mean motion resonances, which are not detected 



in this system ( lAnglada-Escude et al.l . 120121 ) . This black curve is always exterior to the tidal 
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Figure 5: Comparison of the IHZ to the tidal greenhouse limits for Gl 667C c, in a similar format 
as Fig. [2j Here the thin lines correspond to a 4.5 M Earth planet, and thick to 9 M Earth . The 
vertical black lines correspond to the uncertainty in eccentricity of c and d (the latter's existence 
remains uncertain). The dashed, black curve represents where c's and b's orbits cross, and hence 
the region to the left is dynamically unstable. 
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greenhouse curves out to the current orbit of c. Therefore, we may safely conclude that c was 
never in a tidal greenhoi35e_due_tojts orbit . Furthermore, the planet should have become tidally- 
locked within 10 7 years ( IHeller et al.l . 1201 ll ). which is much less than td es - Hence any initial burst 
of tidal heating due to non-equilibrium rotation or obliquity was too short to sterilize this planet. 
Gl 667C c remains a habitable planet candidate. 



7 Discussion 



The previous results show where tidal heating can limit habitability for planets orbiting low mass 
stars. However, many phenomena complicate the process and could alter our findings. Here we 
discuss the results of the previous sections in light of our simplifying assumptions, observational 
requirements, and theoretical limitations. 

We find that Gl 667C c probably did not lose its water to tidal heating because interactions 
with other planets prevent its eccentricity from being large enough to trigger the tidal greenhouse. 
As more potentially habitable planets are discovered around low mass stars, a similar analysis 
as in § should be undertaken in order to assess the possibility that the planet could in fact 
be dehydrated. As we may onl y be able to spectroscopically characterize a few p lanets with the 
James Webb Space Telescope ( ISeager et al.l . 120091 ; Kaltenegger and Traubl . 120091 ) . prioritization 
of targets is crucial, and past and present tidal heating will help determine the best planet to 
observe. 

The timescale for e to decay may be smaller than the timescale for in s i tu terrestrial planet for- 
mation (ILecar and Aarseth . Il986t IWetherill and Stewart! . Il989t iLissauerl . Il993t iKokubo and Ida , 
19981 ; iRaymond et al.l . 120071 ) . On the other hand, terrestri al planets could be pushed into such a 



posit ion by a migrating gas giant in about that timescale ( IRaymond et al.l . l2006t iMandell et al. 



20071 ). It is therefore natural to wonder if such planets can even form, as the tidal effects could 
suppress planet formation, or damp out the eccentricity of a protoplanet before it is massive 
enough to support clement conditions. Several scenarios suggest that fully-formed planets can 
be Tidal Venuses. Orbital instabilities can excite e. Planet-planet scattering and divergent res- 
onan ce crossing appear to play a role in sculpting many planetary systems, in cluding our own 
(e.g. IWeidenschilling and Marzaril . Il996l ; iTsiganis et al.l . 120051 ; iNesvornyl . l201l[ ) . These phenom- 
ena can develop long after planet formation has occurred. Recently it has been suggested that 
exoplanetary systems in resonance appear systematically younger than the general population , 
sugge sting that instabilities can even occur after many billions of years ( IKoriski and Zuckerl . 
20111 ). Therefore, we conclude that tidal greenhouses are plausible. 

The role of oceans in the tidal dissipation process is clearly very important, yet also very 
poorly understood. We have used up-to-date information regarding the dissipation of energy in 
the Earth's ocean and assumed that terrestrial exoplanets will behave similarly, yet many issues 
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remain outstanding. We have ignore d the role of inertia ! waves, which provide additional heating 



of aqueous mantles of icy satellites (jTylerl . 12008k l201lh . Therefore the tidal heating values we 



obtained may, in fact, be too low, further increasing the threat to habitability. Regardless, our 
ranges for Q or r spanned an order of magnitude, but a wider range of values is still possible, 
as they are complex functions of ocean depth, ocean floor topography and the shape of any 
continents or islands that may be present. As "exobathymetry" seems a distant dream, tidal 
dissipation in Earth-like worlds will remain mysterious for the foreseeable future. 

Further complicating the situation, oceans may evaporate long before td es , potentially leading 
to a cyclical process of evaporation and precipitation: After the oceans disappear, Q increases 
and the tidal heat decreases, causing the planet to drop out of the tidal greenhouse, so the 
water rains out, reforming the oceans and lowering Q again. Whether such a cycle exists is 
pure speculation, but we note that an analogous situation can occur with the classic IHZ, where 
evaporated water forms clouds that increase the albedo, which in turn lowers the upward long- 
wavelength flux from the surface, and the planet then drops out of the runaway greenhouse. Our 
choices for Q, r, F cr u and td es can be revised as new observations provide firmer constraints. Of 
these, tdes needs the most work, as it is a fun ction o f poo rly constrained host properties, and 



challenging planetary escape processes, see e.g. iTianl (120091 ) 



A Tidal Venus is an extreme case of tidal heating, over two orders of magnitude more powerful 
than on Io. We have therefore made a considerable extrapolation. Perhaps oceans and/or 
mantles adjust to the increased heating and fail to reach that flux. Our choice for Q an d r imply 



most dissipation occurs in the ocean. For the solid interior, iBehounkova et al.l (120 111 ) find that 
tidal heating at about F crit leads to a "thermal runaway" for planets that transport internal 
energy through convection. Although they did not consider the possibility of advection via 
volcanism, the geophysics of Tidal Venuses require closer scrutiny. We are unaware of any research 
exploring the physical oceanography of planets undergoing strong tidal heating, suggesting it is 
an interesting topic for future research. 

Our treatment ignored mutual gravitational interactions between planets or large satellites, 
which can significantly alter the evolution. Gravitational perturbations can pump up eccentrici- 
ties and obliquities to non-zero values and may be able to modify the spin period. From Fig. [2j 
we see that planets near the inner edge of the HZ may be in a tidal greenhouse for eccentricities 
less than 0.02. Other bodies in the system can easily perturb eccentricities to larger values, 
hence planets in multiple systems may be especially susceptible to a tidal greenhouse, as this 
driven eccentricity can be maintained for arbitrarily long timescales. As individual systems are 
discovered, this point should be revisited. 

On the other hand, orbital stability arguments could allow us to preclude the tidal greenhouse 
in multiplanet systems, as for Gl 667C. When the eccentricity required for the tidal greenhouse 
is so large that the system is unstable, then one can safely exclude the Tidal Venus state. Hence, 
the presence of additional companions provides critical information when assessing habitability. 
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Many other critical phenomena were als o left out such as atmospheric erosion by flaring, stellar 
activ i ty and magnetic dynamo generatio n ( Lammer et al. L 20071; Khodachenko et all 120071 ; iTianl . 
20091 ; ISegura et al.l . l2010t iLammer et all l2010t iDriscoll and Olsonl . 1201 1[ ) . Tidal heating provides 
an interesting counterbalance to atmospheric erosion as it may increase the outgassing rates and 
maintain a permanent atmosphere. The outgassing and escape need to remain in a balance as 
wild pressure and density fluctuations will undoubtedly alter the biosphere, but in principle, 
tidally-driven outgassing could reduce the danger of atmospheric removal. On the other hand, 
more intense outgassing without commensurate draw down by processes like the carbonate- 
silicate cycle could increase the threat of a moist greenhouse. Magnetic fields may slow down 
atmospheric loss, but tidal heating may decrease the dynamo. If magneti c fields are generated by 
convection in the outer core bet ween a hot inner core and a cool mantle (jOlson and Christensenl . 



20061 ; IDriscoll and Olsonl. 1201 11) . then tidal heating of the mantle may suppress magnetic field 



generation ( Stevenson! . 120101 ) . This issue has not yet been explored for tidally- heated exoplanets 



but it could play a major role in habitability. Future work should couple outgassing and escape 
rates in order to determine how the two interact. 



8 Conclusions 



We have shown that tidal heating of some exoplanets may exceed the threshold of the runaway 
greenhouse, the traditional inner edge of the IHZ. We find that for stars with masses <0.3 Ms U n, 
planets in their IHZs with low eccentricity, can be uninhabitable regardless of insolation. We have 
thus fundamentally revised the HZ boundaries for planets on eccentric orbits. Unlike insolation 
from main sequence stars, tidal heating at a desiccating greenhouse level may drop off rapidly, 
but not so rapidly as to preclude the possibility that a planet's entire inventory of water can 
be lost permanently, see Fig. |3j These planets will be uninhabitable regardless of future tidal 
heating, i. e. a planet found with minimal tidal heating today may still have experienced sufficient 
heating for sufficient duration to render it uninhabitable. Additional planetary companions are 
important: They can drive eccentricity and sustain a tidal greenhouse, or they can be used with 
stability arguments to rule out an early tidal greenhouse. 

Traditionally, habitability models have focused on insolation, implying the star is the most 
important aspect of habitability. We have shown that in some circumstances tidal effects are 
more important in determining the inner edge of the HZ. Planetary habit abilit y is a function 
of the star, the planet, and the planetary system (IMeadows et all lin prep.l ). The Kasting et al. 
(119931 ) IHZ has served well as a guide, but is insufficient. Combining all the processes relevant to 
habitability into a single model is a daunting challenge to say the least, but a proper assessment 
of a planet's potential for habitability relies on a wide diversity of properties, some of which will 
not be observable any time soon. Nevertheless, the prospect of identifying an inhabited planet 
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is strong motivation. Moreover, the high cost in time, money and resources required to establish 
a planet as inhabited demand that we use these resources efficiently. In this study we have 
compiled numerous tidal processes and empirical relationships so that at least the tidal effects 
predicted by linear theory may be applied to terrestrial planets in the IHZs of low luminosity 
hosts. 
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A Relations for the inner edge of the habitable zone 

The critical flux for the runaway greenhouse depends weakly on the acceleration due to gravity, 
g, in the planetary atmosphere. Flux is emitted to space from the level of optical depth unity. 
The optical depth of a layer of absorbing gas depends to first approximation on the mass of 
the layer, so with higher g the pressure at the base of a constant mass layer will be higher. 
For the Simpson-Nakajima limit, the T-p structure of the atmosphere is the saturation vapor 
pressure curve, so higher p implies higher T and permits a higher critical flux. This positive 
relationship is offset to some extent by pressure-broadening of absorption, so a t higher q (and 



higher p), less mass of water is re quired to provide an optical depth of unity (jPierrehumbert 
2010l: iGoldblatt and Watsonl I2012J) . 



Pierrehumbertl ( 120101 ) derives a semi-analytical relationship for the critical flux due to the 
Simpson-Nakajima limit: 

/ V 
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where / is the latent heat capacity of water, R is the universal gas constant, P is the pressure 
at which the line strengths are evaluated, g is the gravitational acceleration i n the planetary 



atmos phere, k is a gray absorption coefficient and A is a constant of order unity; iPierrehumbert 
( l2010l ) fits these to numerical runs of a radiative transfer code with n = 0.055 and A = 0.7344 
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where T re f = 273.13 K and P re f = 610.616 Pa are points on the saturation vapor pressure curve 

of water, taken as the triple point. This calculated limit is shown in Fig. [6J 

For convenienc e, we use an empirical relationship for the moist greenhouse limit derived by 

"* (120071 ): 

, lj % \ 1/2 

Hn — \HnQ Oj n i * Oi n l # ) I — I . [o) 

The outer edge (not due to a desiccating greenhouse) can be expressed in a similar format as 

U \ 1/2 

In these equations k n and l out are the inner and outer edges of the IHZ, respectively, in AU, 
U n& and l ou t& are the inner and outer edges of the IHZ in the solar system, respectively, in 
AU, a in = 2.7619 x 1(T 5 AU/K, b in = 3.8095 x 10" 9 AU/K 2 , a out = 1.3786 x 10" 4 AU/K, and 
bout = 1.4286 x 10 -9 AU/K 2 are empirically determined constants, and L* and L are the stellar 
and solar luminosity, respectively. T* = T e ff — 5700 K, where T e ff is the "effective temperature" 
of the star 



boutTDi^- 



T, 



eff 
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(5) 



where R* is the stellar radius and a is the Stefan-Boltzmann constant. Usually, however, the 
IHZ is couched in terms of the stellar mass M* and a, and hence relationships between L*, M* 
and .R* are needed to produce such an expression. Empirical and theoretical relations between 
these properties are reviewed in App. O 

The values of li n& and l ou t& are therefore the key parameters i n the i dentification of the edges 



of the IHZ. We consider three criteria iden tified in ISelsis et all ( 120071 ): 1) 0% cloud cover, 2) 
50% cloud cover, and 3) 100% cloud cover. ISelsis et all ( 120071 ) give values o f l in( ^ of 0.89, 0.72, 



and ~ 0.49 AU for the three possibilities, respectively. For the outer edge, ISelsis et all (120071 ) 



give values of l ou t& of 1.67, 1.95, and 2.4 AU, for the three cloud cover models, respectively. We 
arbitrarily choose the 50% cloud cover case to be the limits of the IHZ. This choice is the middle 
of the possibilit ies, and roughly co rresponds to the cloud cover on the Earth. 



However, as ISelsis et all ( 120071 ) note, "the effect of the spectral type on the albedo, included 
in [these equations] as a quadratic function of (T eff — 5700), was estimated only for a cloud- free 
atmosphere. Since the reflectivity of clouds is less sensitive to wavelength, this quadratic term 
may not be valid to s cale the bou n darie s of the HZ for planets covered by clouds" . Furthermore, 



the assumption that ISelsis et all (120071 ) make that clouds do not affect the infrared flux from 
a planet at the inner edge of the habitable zone is not strong, so the cloudy result s should be 



i nterp reted with caution. Furthermore, assumptions used in the models from which ISelsis et al. 



( 120071 ) derive their results assume fast planetary rotation, so it will not be a good approximation 
for most tidally-locked exoplanets. 
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Figure 6: Surface flux required to trigger a runaway gr eenhouse (F rrii ) as a f unction of planet 
mass. The sol id line is the rela tion for wet planets from iPierrehumbertl (1201 ll ). and the cross is 
the limit from lAbe et al.l (120111 ) for a dry planet. 
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Fig. [7] shows the range of the IHZ for the M dwarf mass range for the empirical relations 
presented in App. For very low-mass stars (< 0.1M Sun ), the position of the IHZ may differ 
by 50% depending on the relationships invoked. However, for larger stars, the IHZ is mostly 
independent of the chosen mass-luminosity and mass-radius relationships. 

B The Desiccation Timescale 



As discussed in § [2], the details of mass loss on a planet near the inner edge of the HZ of a 
very low mass star are largely unknown. Therefore the timescale for desiccation is difficult to 
constrain. In this appendix we describe a simple model for td es and justify the value used in the 
main text of 10 8 years. 



We use the atm ospheric mass l oss mo del described in lErkaev et al.l (J2007I ). which improves the 
standard model by lWatson et al.l (I198ll ). In the Watson et al. picture, a layer in the atmosphere 
exists where absorption of high energy photons heats the particles to a temperature that permits 
Jeans escape. Photons in the X-Ray and extreme ultraviolet have the energy to drive this escape 
on most self-cons i stent atmospheres. In essence, the particles carry away the excess solar energy. 



Watson et al.l (Il98ll ) calculated that Venus would lose its water in 280 Myr. The actual value 



is probably less than that, as they were unaware that XUV emission is larger for younger stars, 
as the spacecraft capable of such observations had yet to be launched. Furthermore, their model 
did not consider th e possibility of mass- loss through L agrange points, or " Roche lobe overflow." 
Erkaev et al.l ( 120071 ) provided a slight modification the IWatson et al.l ( )198ll ) model that accounts 
for this phenomenon for hot Jupiters. The physics of mass loss from a terrestrial atmosphere 
may be different than t hat on a gian t , but as we are interested in objects that are experiencing 
tidal forces, we use the lErkaev et al.l (120071 ) model, 



dM p 



7iR x R p eF X uv 

GM„k tide 



(6) 



where R x is the radius of the atmosphere at which the optical depth for stellar XUV photons is 
unity, e is the efficiency of converting these photons into the kinetic energy of escaping particles, 
Fxuv is the incident flux of the photons. The parameter kude represents the reduction in energy 
required for escape of a particle due to the star's gravitational pull: 
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Figure 7: IHZ boundaries for different combinations of mass- radius and mass-lum inosity relations, 
and using the three different cloud cover IHZ limits from ISelsis et al.l (J20071 ). Except for the 
inner edge near M* = 0.2Ms un , the effective temperature does not affect the boundaries (note 
the numerous curves that are visible). Six combinations are plotted, but most are invisible 
as they are on top of each other. At low M*, the more interior IHZ limits assume the RHOO 
mass-luminosity relationship, see App. O 
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is the ratio of the "Hill radius" to the radius at the absorbing layer. From Eq. 

show that 

Gm p m H k tide 



t, 



ties 



lR„Rn€.F\ 



it is trivial to 



(9) 



where mg is the total mass of all the hydrogen atoms that must be lost. 

The values of Fxuv an d e are therefore the key parameters as they represent the magnitude 
and efficiency of the process. Unfortunately very little is known about the former for very low- 
mass M dwarfs with ages of a few hundred Myr. The stars are very faint, and moreover are 



prone to flaring. However, there is a secular decrease in activity as the stars age (IWest et al 



20081 ). but the data are still too sparse to allow a satisfactory fit. Recent work connected X-Ray 
sources from space telescopes {e.g. GALEX), and ground-based surveys (e.g. Sloan Digital Sky 
Survey), and found tha t older low-mass stars tend to emit about 10~ 6 -10 -3 of their energy in 
the XUV in quiescen ce (Jones and Westl. 120121) . During flares, that emis sion can increase by an 
order of magnitude (JHawley and Pettersenl . Il99ll ; iKowalski et all l2012l ). As all stars are more 
active when they are young, we should expect the XUV emission to be substantially larger at 
those times. 

The value of e must be less than 1, as not all the absorbed energy remov e s par t icles. Observa- 



tions of hot Jupiters suggest val ues of e of .4 f|Yelld . 12004 lLammer et all 120091 ; iJackson et al. 



201fj), while on Venus e ~ 0.15 (jChassefierel . Il996l ). However, our model also requires these pho- 
tons to dissociate the water molecules, hence we should expect e on a Tidal Venus to be much 
less than on a hot Jupiter with a pr edominantly hydrogen atmosphere. For reference, using the 
assumptions of IWatson et al.l (Il98ll ). we find e = 1.7 x 10 -4 implies td es = 280 Myr for Venus. 

In Fig. [8] we show the desiccation timescale for an Earth-like planet orbiting a 0.1 Ms un star 
at 0.03 AU for a range of XUV luminosities and e. The initial mass of the plane t in hydrogen 
i s 1.4 x 10~ 5 M Earth , the current value of the Earth and the same value used in IWatson et al. 
(I198ll ). and R x = 1.02R p . The luminosity of the star is 0.0008 Lg un , and the planet lies inside 
the 0% cloud cover HZ (see App. [A}. The contour lines show logio(tdes)- F° r e J^10~ 4 ; a f ew 
percent of the luminosity would need to be in the XUV for t des to be 10 8 years, while hot Jupiter 
escape rates would only require >10~ 5 of the luminosity to be in the XUV. 

The environment near a forming star is very energetic and therefore we should expect other 
phenomena, such as ablation and magnetic fields, to also erode the atmosphere. Therefore the 
actual values of td es are assuredly shorter than those presented in Fig. [HJ We conclude that 
tdes — 100 Myr is a reasonable choice, while explicitly acknowledging a large range is possible. 

Finally, we note that simulations should couple the evolution of the XU V output to the tida l 
orbital evolution of the planet in order to calculate t^ fo r a given system. IJackson et al.l (120101 ) 



created such a model for CoRoT-7 b fjLeger et al.l . |2009j), which orbits close to solar-mass star. 



For such stars the generic evolution of their XUV evolution is known (jRibas et al.l . 120051 ) and 
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Figure 8: Desiccation timescale for an Earth-like planet orbiting a 0.1 Ms un star as a function 
of the fraction of the luminosity in the XUV and the efficiency of converting that energy into 
escaping particles. Contour lines are logio(td es /yr). The planet initially has the same water mass 
fraction as the Earth and a semi-major axis of 0.03 AU. 
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hence such an analysis was quite revealing. Given the poor constraints on XUV evolution for M 
dwarfs, we eschewed such coupling, deferring that analysis to future work. 

C Masses, Radii and Luminosities of M Dwarfs 

In this appendix we review relationships between mass, radius, luminosity and effective temper- 
ature for M dwarfs (M* < 0.6 Ms un ), and their effects on the limits of the IHZ. Measurements of 
these fundamental stellar properties are challenging and several studies have produced empirical 
and theoretical relations. 



W e first con sider the mass-rad i us rel ation. We include empirical models by lGorda and Svechnikov 
( ll999f ) fGS99). lReid and Hawlevl fepOOh (R HOOt. lBavless and Oroszl fl2006h rBO06^ as well as the 
theoretical models of iBaraffe et al.l ( 19981 ) (B98). This latter model is parametrized in terms of 
the metallicity as defined by astronomers: 



[Fe/H] = log 10 [Fe/H], - log 10 [Fe/H] 



Suni 



(10) 



in other words it is the difference between the ratio of iron to hydrogen in a star to that of the 
sun. B98 consider solar metallicity stars ([Fe/H] = 0) and sub-solar ([Fe/H] = -0.5). GS99 and 
BO06 present analytic fits that are shown in Table 1. For RH00, we take their data in Table 4.1 
and derive a third-order polynomial fit, y = b + b\X + b 2 x 2 + b 3 x 3 , using a Levenberg-Marquardt 
minimization scheme, and present the results in Table 1. The top two panels of Fig. |H]show the 
relationships graphically. 

Next we examine the mass-luminosity relationships detailed in RHOO. IScalo et al.l (120071 ) (S07) 
and B98. For the former, we again fit their data with a third-order polynomial. The fits are 
listed in Table 1, and shown graphically in the bottom panels of Fig. ED 
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Figure 9: Scaling relations for M dwarfs. Top: Mass-radius relations. Black curves are empirical 
relations: solid from BO06, dashed from RHOO, and dotted from GS99. Red curves are theoretical 
curves from B98: solid assumes [Fe/H] = 0, dashed [Fe/H] = -0.5. Bottom: Mas s-luminosity 



relatio ns. Black curves are empirical relations: solid from RHOO, and dotted from IScalo et al. 



( 120071 ). Red curves are theoretical curves from B98: solid assumes [Fe/H] = 0, dashed [Fe/H] 
-0.5. 
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Tabic 1: Scaling Relations for M Dwarfs 
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D Radii of Terrestrial Planets 



In this appendix we review published studies which provide analytic formulae relating terrestrial 
planet mas s and radius. Throughout the study we appealed to the relationship derived by 
Sotin et al.l ( 120071 ). who considered two types of planets: Earth-like composition and "ocean 
planets" that are 50% H2O by weight. The scaling relationship for the former is 
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(11) 



(12) 



Fortney et al.l ( 120071 ) parametrize the radii (M p > 0.01 M Ea rth) in terms of the ratio of ice to 
rock, f ice , or rock to iron, f roc k- They assume that these two combinations are the most likely 
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for terrestrial exoplanets and derived the following relationships: 
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(0.0912/ Jce +0.1603)(log 10 J^) 2 +(0.330/ ice +0.7387) log 10 -^— 

JVlEaxth MEarth 



(0.0592/ rocfc +0.0975)(log 10 J^) 2 +(0.2337/ rocfc +0.4938) log 10 ~^- 

JViEarth MEarth 



f-0.4639/ ice +1.1193, 
(13) 
f-0.3102/ rocfc +0.7932. 

(14) 



In these equation f ice = 1 corresponds to a pure iceball, fi ce = and f roc t = 1 to a pure rock 
(silicate) planet, and f roc k = is a planet made of pure iron. Note that Eqs. (fT3"l) - (fT4"l) are 
analytic fits to complex models, and hence predict slightly different values for R p for a planet 
made of rock (f ice = and f roc k = !)• The Earth has f roc k = 0.67 and hence Eq. ( TT41) becomes 
R v l Rfarth = 0.1372[login (Mp/ M EaJ th)] 2 + 0.6504 log 10 (M p /M Ea rth) + 1.0010. 

Valencia et al.l (120071 ) consi dered planets i n the range 1-10 MEarth and parametrized the 
radius in terms of fi ce (see also IValencia et al.l . 120101 ). They find 



Ri, 



R 



Earth 



;i + 0.56/,, 



M n 



0.262(1-0. 138/ lce ) 



M 



Earth 



(15) 



Seager et al.l (120071 ) consider a wide range of planetary compositions and develop a general, 



but complicated, mass-radius relationship. Although they developed several formulae, here we 
only present one. They found that by casting M p and R p in terms of "scaling values" all planets 
follow the following relationship: 



log 
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(16) 



where k\ = —0.20945, k% = 0.0804, and k% = 0.394. The parameters m s and r s are the scaled mass 
and radii, respectively, and are defined as m s = m/rrti and r s = r/ri, where the denominators 
are a scaling factor that depends on composition and thermal effects. We will consider four 
compositions. For pure iron planets, m\ = 5.8 M Ea rth, and r 1 = 2.52 Rsarth- For pure rock 
(perovskite) planets, mi = 10.55 M Ear th, and r\ = 3.9 REarth- For pure water planets, m 1 = 
5.52 M Ear th and r% = 4.43 REarth- Finally, for an Earth-like planet (30% iron, 70% perovskite), 
mi = 6.41 MEarth and r% = 2.84 REarth- No te that this formulat ion is only valid for m s < 4. 

The last relationship we review is that of iGrasset et al.l ( 120091 ). who cast their model in terms 
only of fi ce . Their formulation is complicated and rather than reproduce it here, the r e ader i s 
referred to their paper. We do not include the recent models put forth by ISwift et al.l ( 120121 ). 
as they do not incl ude f ormulae, but they do find their results are consistent with those of 
Fortney et al.l (120071 ) and [Valencia et al.l ( 120101 ) for rock and iron planets. 

In Fig. [10] we compare these relationships from 0.3 - 30 MEarth? and for a variety of compo- 
sitions. Not all the relationships span that mass range, nor do they encompass all compositions. 
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Furthermore, a lthough we group p lanets into four types, not all types are identical. For example, 
we labeled the ISotin et al.l ( 120071 ) ocean planet as a "Water" planet, but it is only 50% water, 
while the other three are 100% water. Nonetheless, the plot shows that if only mass is known, 
the radius of a planet can only be known to a factor of 2-2.5 (of course, its terrestrial nature 
is also unlikely to be known). As tidal heating scales as Rp, this translates to a difference in 
heati ng of 30-100, assu ming identical Q or r. Note that the solid red curve (Earth-like planet 
from lSotin et al.l ( 20071 )) has a typical shape for its class, hence the results presented here do not 
vary much if a different model is chose, unless a significantly different composition is invoked. 



E Tidal Theory 

For our c alculations of tidal evolution, we employ "equilibrium tide" models, originally con- 



ceived by iDarwinl (11880l ). This model assumes the gravitational potential of the tide raiser can 
be expressed as the sum of Legendre polynomials (i.e. surface waves) and that the elongated 
equilibrium shape of the perturbed body is slightly misaligned with respect to the line which 
connects the two centers of mass. This misalignment is due to dissipative processes within the 
deformed body and leads to a secular evolution of the orbit as well as the spin angular momenta 
of the two bodies. As described below, this approach leads to a set of 6 coupled, non-linear 
differential equations, but note that the model is, in fact, linear in the sense that there is no 
coupling between the surface waves which sum to the equilibrium shape . A substantial body 



of research is devoted to tidal t heory (e.g. iHutl. Il981t iFerraz-Mello et al.l . 120081 ; IWisdoml . 12008 



Efroimsky and Williamsl . 120091 ; iLeconte et al.l . |2010| ). and the reader is referred to these studies 



for a more complete description of the derivations an d nuances of tidal t heory. For this inves- 
tigation, we will use the models and nomenclature of iHeller et al.l ( ]201ll ). which are presented 
below. 



E.l The Constant Phase Lag Model 

In the "constant-phase-lag" (CPL) model of tidal evolution, the angle between the line connecting 
the centers of mass and the tidal bulge is assumed to be constant. This approach has the 
advantage of being analogous to a damped driven har monic oscillator, a well-studied system , 
and i s quite commonly utilized in planetary studies (e.g. iGoldreich and Soterl . Il966l ; iGreenbergl . 
20091 ). In this case, the evolution is described by the following equations 
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Figure 10: Mass-radius relationships for terrestrial planets of varying compositions. Solid curves 
assume a roughly Earth-like composition, dashed are all silicates, do t ted ar e at least half water, 

and d ot-dashe d are pure iro n. Black corres pond to iFortnev et al.l (120071 ). red t o ISotin et al. 
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where e is eccentricity, t is time, a is semi-major axis, G is Newton's gravitational constant, Mi 
and M 2 are the two masses, R\ and i? 2 are the two radii, Ui are the rotational frequencies, ipi 
are the obliquities, and n is the mean motion. The quantity Z[ is 
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(19) 



where /c 2) j are the Love numbers of order 2, and Qj are the "tidal quality factors." The parameter 
& is 
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(20) 



where i and j refer to the two bodies, and r g is the "radius of gyration," i.e. the moment of 
inertia is M(r g R) 2 . The signs of the phase lags are 



(21) 



£ ,i = S(2wj - 2n) 
£ 1; j = S(2cjj — 3n) 
£ 2 ,i = S(2wj - n) 
e 5>i = S(n) 
£ 8ii = S(wj - 2n) 
£9,i = S(wj) , 

with E(x) the sign of any physical quantity x, thus S(x) = + 1 V — 1 V 0. 

The tidal heating of the ith body is due to the transformation of rotational and/or orbital 
energy into frictional heating. The heating from the orbit is 
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and that from the rotation is 
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The total heat in the zth body is therefore 
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The rate of evolution and amount of heating are set by three free parameters: Q,k2, and r g 
dis cuss the ranges o f these parameters below in App. IE14. 



Goldreichl (119661 ) suggested that the equilibrium rotation period for both bodies is 



dCPL 
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Murray and Dermottl (119991 ) present a derivation of this expression, which assumes the rotation 
rate may take a continuum of values. However, the CPL model described above only permits 
4 "tidal waves" , and hence does not permit this continuum. In Fig. [TTJ we compare ( 1231) to the 
equilibrium values predicted by Eqs. (ITTj)-(l2~T]) as a function of eccentricity for Gl 581 d , i.e. an 



orbital period of 66 days, a mass of 5.6 Msarth and a radius of 1.6 R.Earth (jMayor et al.l . 120091 ). 
As expected (123]) predicts a continuous range of periods, whereas our model (solid curve) does 
not. The jump at e = yl/19 occurs at the ui/n =1.5 transition, and nearly reaches the value 
predicted by (123]) . The next phase jump at e = a/2/19 ~ 0.32 is not included in the CPL model, 
and hence it breaks down at larger values. However, at larger e values, coupling between waves 
and other non-linear effects may dampen the role of tides. Therefore the evolution at larger e 
predicted by the CPL model may not be qualitatively correct. We urge caution when interpreting 
CPL results above e = 0.32. 

Thus, the CPL model can be considered from two perspectives. When modeling the evolu- 
tion of a system, one should use the discrete spin values for self-consistency, i.e. as an initial 
condition, or if forcing the spin to remain tide-locked. However, if calculating the equilibrium 
spin period separately, the continuous value of Eq. [23] should be used. We refer to these rotations 
as "discrete states" and the "continuous state" . Note that thi s point mark s a diff erence between 
our implementation of the CPL model, and that described in lHeller et al.l (J20111 ). 



E.2 The Constant Time Lag Model 

The constant-time-lag (CTL) model assumes that the time interval between the passage of the 
perturber and the tidal bulge is constant. This assumption allows the tidal response to be 
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Figure 11: Comparison of the equilibrium spin periods in the CPL framework. The solid line 
is the value pre dicted by our model, Eqs. (TT71)-(12T1). the dashed is the "continuous" model of 
Goldreichl ( fl966h . 
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continuous over a wide range of frequencies, unlike the CPL model. But, if the phase lag is a 
function of the forcing frequency, then the system is no longer analogous to a damped driven 
har monic oscil l ator. T herefore, this model should only be used over a narrow range of frequencies, 
see iGreenbergi (120091 ) . However, this model's use is widespread, especially at high e, so we use 



it to evaluate tidal effects as well. This model predicts larger tidal heating and evolution rates 
at high e as any coupling between tidal waves is ignored. Therefore, the CPL and CTL models 
probably bracket the actual evolution. 

The evolution is described by the following equations: 
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The tidal heating of the ith body is therefore 
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Note that Eqs. (I2"2"j) and (1501) include terms due to both rotation and obliquity, which were not 
discussed in the main text. Should information regarding these parameters become available for 
particular planets, they should be included in estimates of tidal heating and habitability. 
It can also be shown that the equilibrium rotation period for both bodies is 



q 2/ 2 (e) cos ip 

which for low e and tp = reduces to 

PS,™ = T^r (32) 

There is no general conversion between Q p and r p . Only if e = and ip p = 0, when merely 
a single tidal lag angle e p exists, then 

Q p « l/(2|n-u;p|Tp), (33) 

as long as n - Up remains unchanged. Hence, a dissipation value for an Earth-like planet of 
Q p = 100 is not necessarily equivalent to a tidal time lag of 638 s, so the results for the tidal 
evolution will intrinsically differ among the CPL and the CTL model. However, both choices are 
common for the respective model. 

E.3 Numerical Methods 

The rates of change of the 6 parameters affected by tidal theory are all different. The orbital 
parameters tend to evolve several orders of magnitude more slowly than the planetary spin 
parameters, at least for the cases we are interested in. Therefore we implement a dynamical 
timestep scheme in our integration of tidal evolution. At the beginning of each timestep, we 
calculate the derivatives and set the timestep to 

At = 77 x minf -?—, -?—, U * A Up ^ ) , (34) 

\da/dt' de/dt : duj*/dV dip^/dt' du p /dt : dip p /dt/' 

where 77 is a constant, which we set to 0.01. We find that larger values of 77, even up to unity, 
produce results similar with 77 = 0.01, but at our more conservative choice, we find almost no 
difference with even smaller value of rj, i.e. the solution has converged. This implementation 
increases the speed of the code by several orders of magnitude in some cases. For example, if 
considering a planet with spin properties similar to the Earth at 0.01 AU, the timescale to damp 
out the obliquity is of order years, but the orbit evolves in millions of years (depending on e). 
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Integrating this system for billions of years with a fixed timestep chosen to resolve the obliquity 
evolution will take many orders of magnitude more time than using this sche me. Dynamical 



timestepping is especially important for planets that merge with their host star ([Jackson et al 



20091 : iLevrard et all 120091 ) 



In many cases, e and ip can damp to arbitrarily small values, and hence de/dt and dip/dt take 
on meaningless values. In some cases, their respective timescales as defined above can become 
very short, even though they are effectively 0. This situation destroys the speed advantage gained 
by utilizing (IM|) . We therefore set a floor for these values at 10~ 10 . When they reach this level, 
we set them to 0, and do not include them in the determination of At. 

Finally, care must be taken as a planet approaches tidal locking. In many cases an integration 
will bounce back and forth above and below the equilibrium period. This possibility is a natural 
result of the discrete nature of numerical integration, but, again, can diminish the benefits of a 
dynamical timestep. Therefore, if a planet comes within 10% of the equilibrium period, we set 
the spin period to the equilibrium value. If the orbit is eccentric and evolving, after tidal locking, 
we force the planet to remain locked, i.e. the spin period obeys Eqs. ( 123]) or (l3Tj) . 



E.4 Tidal Response in Celestial Bodies 

Exoplanets are expected to hav e an even wider diversity of composition s and structure s than 



those in our Sola r System (see iRaymond et all I2004J ; iLeger et all 12004 ; lO'Brien et all . 12006 



Bond et all l2010l ). thus their response to tidal processes probably varies greatly. We therefore 
choose to keep our model simple and set k 2 = r g = 0.5 for all planets. 

The Earth is the only body in the Solar System with a significant amount of water on its 
surface, and is also the planet with the most detailed picture of its tidal processes. The majority 
of tidal dissipation on the Earth occurs in the oceans, and we assume the same will occur on 
terrestrial exoplanets. Tidal dis sipation in the Ea r th's o cean was directly measured from TOPEX 
and Poseiden satellite data by Egbert and Ray! (120001 ) who found that 50-75% of dissipation 
occurs in and around shallow seas, with the remainder in the open ocean, presumably due to 
turbulence c aused by topography. Recent computer simulation s have successfully reproduced this 
observation (jJaynd . l200ll; iNvcanderl. 12 005; K elly et all 120101 ) . For a more complete treatment 
of this process, consult [Garrett and Kunzei ( 120071 ). Therefore dissipation in planetary oceans is 
a complex function of continental shapes, sea floor topography and the forcing frequency. 

Lunar Laser Ranging, in which the timing of laser pulses bounced off the Moon is measured, 
has provided very p recise measurements of the rate of the Moon's recession from the Earth 



( IDickey et all Il994f ). If this recession is d ue only to tid al dissipation within the Earth, than the 
Earth's tidal quality factor is Q Earth = 12 (JYoderl . Il995l ). This value is lower than most measure- 
ments in the Solar System, which are more uncertain, see below. However, given the presence 
of surface oceans, this extra dissipation should be expected. On the other hand, extrapolating 
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backward from the current orbital expansion ra te and employing t he CPL model, one finds the 
Moon was at the Earth's surface ~ 2 Gyr ago (JMacDonaldl . Il964l ). The moon-forming impact 
has been dated to about 4.5 Gyr ago, and hence the Earth's current Q appears too low. This 
discrepancy may be explained by 4 possibilities: 1) The current Q is anomalously low, and in 
the past different continental configurations and/or ocean floor topography did not permit such 
a high dissipation, 2) The tidal effects are a function of frequency, and as the Earth's rotation 
has slowed , the dissip ation rate has increased, 3) Perturbations from other planets affected the 
evolution ( iCukl . 120071 ). or 4) Linear theory is oversimplified. Probably all of these possibilities 
play a role, and make the extension of tidal theory to an Earth-like exoplanet unreli able. 

Measurements of the Martia n dissipation functi on find QMars = 85.58 ± 0.37 (IBills et al. 
20051 ) or QMars = 79.91 ± 0.69 (ILainev et all 120071) . The value for Mercury is (^Mercury < 190, 
and for Venus it is Qy en „ s < 17 (IGoldreich and Soterl . Il966[ ) . Observations of Io im ply Qto < 1 00 



( Peale et all Il979l ; ISegatz et all . Il988h . while the Moon lies at Q 



/Moon 



26.5 ±1 flDickev et al 



19941 ) . These measurements indicate a similar order of magnitude for the tidal dissipation func- 
tion of all desiccated bodies in the Solar System. 

Tidal effects on Venus may be particularly complicated due to its large atmosphere. The 
rotation rate of Venus is such that it always presents the same face to the Earth at inferior 
conjunction. Although this may be by chance due to Venus approaching tidal lock with the Sun, 
it has also been suggested that at some point in the past thermal atmospheric tides could have 
provided an equal but opposite torque than that from the rocky component. If that o ccurred, 
the torque from the Earth could dominate and spin-lock Venus (IGold and Soterl . I19691 ). More 
recent and detailed work has shown that this interaction can be sig nificant and permits chaoti c 
spin evolution that eventually settles into 1 of 4 possible states ( ICorreia and Laskarl . 120011 ). 
Although a large atmosph ere can play a significant role in the evolution of terrestrial planets 
( ICorreia and Laskarl . 120031 ). we ignore their effect here due to the large number of unknowns for 
an exoplanet. 

In real bodies, Q is a function of the b odies' rigidity //, viscosity 77, and temperature T 
( ISegatz et all Il988t iFischer and Spohnl . Il990l ) . A comprehensive tidal model would couple the or- 
bital and st ructural evolution of the bodies since small perturbations in T can result in large varia- 
tions in Q (ISegatz et all Il988l ; iMardling and Linl . |2002| ; lEfroimsky and Laineyl . 120071 ) . Moreover, 
tidally-generated heat will not necessarily be transported to the surface quickly. The internal 
structure of the planet may trap the energy for long periods of time, resultin g in catas t rophi c 



tectonics when the heat is finally released, as may have occurred on Venus (jTurcottel . 119961 ). 
This episodic volcanism could change td es dramatically, but, on the other hand, complete litho- 
spheric overturn is probably j ust as detrimen t al to life as desiccation. Nonetheless, future work 
should explore these details. iHenning et al.l ( 120091 ) performed detailed calculations and argue 
that k 2 = 0.3 and Q = 50 are reasonable choices for dry planets. 

The tidal properties of stars are about as poorly constrained as for planets. The primary 
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constraint comes from the distribution of orbits of binary stars (e.g. IZahnl . 119771 ). whose orbits 
are circularized by the tidal interaction. However, for the closest binaries, where the tidal effects 
are strongest, and one might think th e constraints the tightest , the situation is complicated by 
the early phases of stellar evolution ( IZahn and Bouchetl . Il989l ). in which stellar radii contract 
rapidly As the radius e nters the tidal evolution a t the fifth power, the radial contraction is 
a major effect. Recently iKhaliullin and Khaliullinal (J201ll ) revisited this point and found that, 
although radial contraction does play a major role, they were unable to comb i ne mo dels of the 
equilibrium tide and radial contraction to match observations. IWeinberg et al.l (j201l[ ) tested the 
linear theory against non-linear models and found that in general the linear theory does not 
capture the physics of the more sophisticated approach. These studies demonstrate that there 
still exist serious problem with the linear theory, and quantitative discrepancies exist. However, 
it does qualitatively model the relevant physics, and allows wide sweeps of parameter space to be 
a computationally tractable problem. We therefore continue to use it, but explicitly acknowledge 
that known weaknesses persist. 

A wide range of values for the t idal quality factor a nd tidal time lag for star s are present in the 
litera ture, ranging from Q* = 10 4 ( lAdams et al.l . 1201 ll ) all t he way toQ, = 10 9 (IMatsumura et al. 



2010) , with many others finding intermediate values (e.g. iLin et al.l . Il996l ; ICarone and Patzoldl . 



20071 ; iJackson et al.l . l2008bl . 120091 ; llbgui and Burrows! . 120091 ) . For the star we take an intermediate 
and typical value of Q* = 10 6 . 
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